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Genome Sequence of the Asiatic Species Borrelia persica 
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We report the complete genome sequence of Borrelia persica, the causative agent of tick-borne relapsing fever borreliosis on the 
Asian continent. Its genome of 1,784,979 bp contains 1,850 open reading frames, three ribosomal RNAs, and 32 tRNAs. One clus- 
tered regularly interspaced short palindromic repeat (CRISPR) was detected. 
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Borrelia persica is the agent of relapsing fever (RF) in central 
Asia and Middle Eastern countries, but lack of awareness 
about the disease probably has led to an underestimation of its 
prevalence (1). Until now, B. persica has eluded in vitro cultiva- 
tion, impeding its full characterization and in vitro experiments. 

The B. persica strain No 14 was isolated from Ornithodoros papil- 
lipes collected near the city of Khiva, Uzbekistan. The strain was 
isolated and propagated at the Laboratory of Vectors of Infections, 
Gamaleya Research Institute, Moscow, Russia, by passaging in 
guinea pigs (2) . It was initially kept alive in a tick-guinea pig cycle and 
in 1988 it was stored at -70°C. After thawing in 2008, the strain was 
refractory to growth in vitro and was therefore propagated in CB17/ 
Icr-Prkdc''"''/IcrIcoCrl mice (Taconic, Denmark) and subsequendy 
successfully grown at 37°C in freshly made Barbour-Stoenner-Kelly 
II (BSK-II) medium with 1.4% (wt/vol) gelatin and 10% (vol/vol) 
rabbit serum (3). DNA was extracted using the Wizard genomic 
DNA purification kit (Promega Biotech AB, Sweden) and sequenced 
using an lUumina HiSeq 2000 sequencer. De novo assembly of reads 
was performed with the ABySS 1.3.4 assembler, and open reading 
frames (ORFs) were identified with Prodigal (4). tRNAs were pre- 
dicted with the Aragorn software (5), and ribosomal RNAs were pre- 
dicted with RNAmmer. The genes were farther annotated using 
BLAST against the NCBI nonredundant database. The functional 
categorization of ORFs was performed using online RPS-BLAST 
against the cluster of orthologous groups (COG) (6) and Pfam data- 
bases (7). PHAST (8) and the Prophage Finder software were used for 
bacteriophage detection (9). Clustered regularly interspaced short 
palindromic repeats (CRISPRs) were predicted using the CRISP- 
Rfinder server (http://crispr.u-psud.fr/Server/). To estimate the sim- 
ilarity with pubHshed RF Borrelia genomes (10, 11), the average nu- 
cleotide identity (ANI) was calculated as described in Konstantinidis 
etal. (12). 

A total of 1,795,421 bp of sequence was generated from 
15,331,580 reads with 826-fold coverage and a G+C content of 
28.7%. The scaffolds were ordered, resulting in a linear chromo- 
some with a size of923,419bp and unclosed plasmids with a size of 
872,002 bp. Gene finding resulted in 1,850 ORFs, three ribosomal 
RNAs, and 32 transfer RNAs, which is similar to the genomes of 



other RF species (10, 11). The protein-coding regions represent 
81% of the B. persica genome. Among the predicted genes, 43.4% 
were assigned a COG function. One CRISPR with three spacers 
was found. Comparing the chromosome of B. persica to the pre- 
viously published genomes of Borrelia duttonii, Borrelia recurren- 
tis, and Borrelia crocidurae revealed a high colinearity, with an 
average nucleotide identity of 89%. 

Several genome sequences of African RF Borrelia spp. have 
been published (10, 11). The added B. persica genome sequence is 
of great importance as the first representative of Asian tick- 
borne RF borreliosis. An analysis of its genome sequence will pro- 
vide further insights into its functional genomics and host adap- 
tation. It wiU also facilitate future molecular and medical research 
on B. persica, a bacterium causing human disease, about which 
very little is known. 

Nucleotide sequence accession numbers. This whole-genome 
shotgun project has been deposited at DDBJ/EMBL/GenBank un- 
der the accession no. AYOTOOOOOOOO. The version described in 
this paper is version AYOTO 1000000. 
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